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Abstract. We demonstrate the in situ detection of cold S7 Rb atoms near a dielectric surface using the 
absorption of a weak, resonant evanescent wave. We have used this technique in time of flight experiments 
determining the density of atoms falling on the surface. A quantitative understanding of the measured 
curve was obtained using a detailed calculation of the evanescent intensity distribution. We have also 
used it to detect atoms trapped near the surface in a standing-wave optical dipole potential. This trap 
was loaded by inelastic bouncing on a strong, repulsive evanescent potential. We estimate that we trap 
1.5 x 10 4 atoms at a density 100 times higher than the falling atoms. 
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PACS. 32.80.Pj Optical cooling of atoms; trapping 
sorption 



1 Introduction 



Recently there has been increased interest in cold atoms 
trapped near a surface. For example magnetic fields of 
micro-electronic structures are used to trap and guide 
atoms near the surface of so-called atom chips [1-3]. Last 
year Bose-Einstein condensation was realized on such a 
chip [4,5]. Other examples are experiments aiming to in- 
crease the phase space density and create a quasi 2D gas of 
atoms using inelastic reflections from an evanescent-wave 
(EW) atomic mirror [6-8]. 

These experiments pose new challenges for in situ de- 
tection, in particular if the atoms are within the order of 
an optical wavelength from the surface. In this case ab- 
sorption of EWs is advantageous, since only atoms that 
are within a few wavelengths from the surface are probed. 
Aspect et al. have proposed the nondestructive detection 
of atoms close to a surface by detecting a phase change 
of a far detuned EW [9]. However this effect is very small 
and has not been observed so far. In this letter we demon- 
strate experimentally the absorption of resonant evanes- 
cent waves as a novel diagnostic tool to study cold atoms 
near a dielectric surface. Using a weak, resonant EW, we 
have studied a sample of cold atoms falling onto the sur- 
face as well as atoms trapped close to the surface. 

EW absorption has previously been used for spectros- 
copy on hot and dense atomic vapors to experimentally 
investigate EW properties such as the Goos-Hanchen shift 
[10,11]. 



42.25.Bs Wave propagation, transmission and ab- 



2 Evanescent wave calculations 

An evanescent wave appears when a light wave undergoes 
total internal reflection at the surface of a dielectric with 
index of refraction n. For a plane incident wave the optical 
field on the vacuum side of the surface decays exponen- 
tially ~ exp(— z/£) with z the direction perpendicular to 
the surface, £(0) = ^(n 2 sin 2 9 - l) -1 / 2 the decay length, 
n the index of refraction of the substrate and 9 the angle 
of incidence. For a low density of resonant absorbers near 
the surface, the scattering rate in the low saturation limit 
is proportional to the square of the field: ~ exp(— 2z/£). 
If the density of absorbers is uniform, this gives rise to 
a total rate of scattered photons proportional to £. The 
scattered photons are missing from the reflected beam, 
which is therefore attenuated. If the angle of incidence ap- 
proaches the critical angle 9 C = arcsin(n _1 ), the value of £ 
diverges, so the absorption becomes large. The absorption 
is however less height selective in this limit. 

For a Gaussian beam with an angle of incidence 9q ^> 
9 C + <t> with 9 C the critical angle and <j> the divergence of 
the beam, the electric field is similar to the field of a plane 
wave. For Gaussian beams with 9o closer to 9 C , the evanes- 
cent field does not decay as a simple exponential. We can 
describe the incident field as a superposition of plane wave 
Fourier components with various 9. Each component con- 
tributes an evanescent field component with decay length 
£(0) and an amplitude proportional to the complex Fresnel 
transmission coefficient t(0). Because both £(0) and t(9) 
vary strongly around 9 C , the evanescent wave contribu- 
tions of these incident plane wave components add up to 
a non exponentially decaying field. In addition the trans- 
verse beam profile is modified to a non Gaussian shape. 
In the reflected beam one of the effects is the well known 
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Fig. 1. (a) Overview of the geometry and notations for the 
evanescent wave calculations near the critical angle. The inci- 
dent beam has a waist w and angle of incidence 6q . The waist 
is at a distance d from the surface, (b) Intensity distribution 
of the evanescent wave for realistic experimental parameters 
w = 330 lira, (6o — 8 C ) = 133 ^irad and d = 680 mm. 
Transverse i-distribution at the prism surface, z-distribution 
at the x-coordinate where the intensity at the surface is high- 
est. 



Goos-Hanchen shift [12]. Other phenomena like nonspecu- 
lar reflection, shift of focus point and change of beam waist 
have been predicted [13]. They all result from combining a 
finite size beam with angle dependence of the transmission 
coefficient t(9). Recently, it has been proposed to measure 
a Goos-Hanchen shift also in the evanescent wave using a 
scanning tunnelling optical microscope (STOM) [14]. 

In the following calculations we will coherently add 
plane wave components with propagation vectors in the 
plane of incidence. The transverse distribution has been 
assumed Gaussian with the same radius as the incident 
beam at the surface. This approach is valid, since the 
transverse components of the propagation vector change 
the angle of incidence only in second order. 

The evanescent field has been calculated by evaluating 
the following expression 



tv/2 



E(x, z) 



ts, P (0)p(0) exp[ik x (9)x 



+l nk4(6-6 y -i£z*£]d0, (1) 



with k x (9) = nko sin 9 the wavevector along the prism sur- 
face, fcrj = 27r/A and 9 is the angle of incidence. The first 
two exponents are the field distributions parallel and per- 
pendicular to the surface, respectively. The third exponent 
takes into account that the waist of the incident beam is 
at a distance d from the surface. The fourth exponent is 
the distribution over angles of incidence for a Gaussian 
beam, with <j> = 2/(nkow) the 1/e half width of the angle 
distribution of the field of a Gaussian beam with waist w. 
The factor t s>p (8) is the Fresnel transmission coefficient for 
transmission of a s or p polarized plane wave with angle 
of incidence 9 from the dielectric into vacuum. They are 
given by t s (9) = 2n cos 9 / (n cos 9 + i\/ n 2 sin 2 9 — 1) and 

t p (9) — 2n cos 9/ (cos 9+in\/ n 2 sin 2 9—1) respectively. Fi- 
nally p(9) is a normalization factor that is equal to 1 for 

s polarized incident light and \J 2n 2 sin 2 9—1 for p polar- 
ized light. The integration is carried out over the range of 
angles of incidence that generate an evanescent wave, from 
the critical angle 9 C to n/2. The normalization is chosen 
such that |.E| 2 = 1 in the waist of the incident beam. The 
geometry of the problem and some of the parameters are 
displayed in Fig. la. The effective evanescent intensity is 
given by 

I(x,y,z) = -I Q — \E(x,z)\ 2 e- 2 y 2 H, (2) 
n w p 

where Iq is the intensity of the incident beam in the waist 
of the incident beam. The Gaussian determines the distri- 
bution in the y-direction with w p the transverse 1/e 2 in- 
tensity radius at the prism surface. The fraction w/w p cor- 
rects for the transverse divergence of the incident beam. 
This approach is possible since \E\ 2 — 1 in the waist of the 
incident beam. Fig. lb shows a calculation of the EW in- 
tensity distribution for realistic experimental parameters 
($o — #c) = 133 /zrad, w = 330 /im and d = 680 mm. 
This distribution is used to calculate the absorption via 
the integrated scattering rate. 



3 Time of flight experiments 

In a first experiment we use evanescent wave absorption 
to detect a cloud of cold atoms falling on a glass surface. 
Our setup consists of a right angle, uncoated, BK7 prism, 
mounted inside a vacuum cell with a background Rb pres- 
sure of 10~ 7 - 10 -8 mbar. About 2 x 10 7 87 Rb atoms 
are captured from the background gas in a magneto op- 
tical trap (MOT), located 7 mm above the prism surface 
and are subsequently cooled to 4.5 fiK in optical molasses. 
They are released in the F=2 hyperfine ground state and 
fall towards the prism due to gravity. Inside the prism a 
weak, resonant, p-polarized probe beam undergoes total 
internal reflection. Its angle of incidence was determined 
with 100 /irad accuracy by a method that will be discussed 
later. The probe beam has a waist of 330 /im (1/e 2 radius) 
resulting in a divergence of 500 /xrad inside the prism. At 
the prism surface it has a waist (1/e 2 intensity radius) of 
(470 ± 20) /xm. The total power of the incident beam is 
P in = 2.2 fJ,W. 
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Fig. 2. (a) measured absorption by a cloud of cold atoms 
falling onto a dielectric surface as a function of the angle of 
incidence of the evanescent probe beam. The inset shows a 
typical time of flight signal, measured with an angle of inci- 
dence of (80 — 9 C ) = 130 fj,ra,d. (b) same data as shown in 
(a), but versus effective decay length £. Also the results of the 
calculations based on a density of 1.2 x 10 9 cm -3 are shown. ■ 
Measured data, O calculated absorption on the basis of scat- 
tering, A calculations using complex index of refraction. The 
vertical dotted line corresponds to the critical angle. 



We measure the absorption in the reflected wave on a 
photodiode which is amplified by a factor 10 6 by a fast, 
low noise current amplifier (Femto HCA-1M-1M). A typ- 
ical absorption trace is shown in the inset of Fig. la. The 
maximum absorption of time traces for different values of 
the angle of incidence are plotted in Fig. 2a. From this 
graph it is clear that the absorption is highest for angles 
of incidence very close to the critical angle. In order to 
analyze these results we consider the atomic density to 
be uniform perpendicular to the surface, since the pene- 
tration of the EW (< 10 /im) is much smaller than the 
size of the cloud (~ 1 mm). It is crucial to take the fi- 
nite probe beam diameter into account. This leads to a 
finite range of angles of incidence so that the EW is no 
longer exponential as described above. We define an ef- 
fective decay length £ by £/2 = J (z 2 ) — (z) 2 where the 
distribution function is the normalized intensity distribu- 
tion I(xq, 0, z)j J Q I(xq, 0, z)dz at the transverse position 
Xq where the intensity at the surface is maximum. For a 
plane incident wave £ = £. In Fig. 2b the solid squares 
represent the same absorption data as shown in Fig. 2a, 
but plotted versus £. Absorption increases with £, but de- 



creases beyond a value £ w 12A. This decrease for larger 
£ occurs because the amplitude of the EW quickly drops 
for angles of incidence 9q smaller than 9 C . 

We compare our data to the absorption as calculated 
using two different approaches. The first method is to cal- 
culate the scattering of evanescent photons by the atoms 
near the surface, where we assume that the transmission 
coefficients are not changed by the atoms. The scattered 
power is calculated as 



1 



-hwr 



S(x) 3 



EW 



a(x) 
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where s(x) = ^/(x)/J sat is the local saturation param- 
eter, 7 sa t is the saturation intensity 1.6 mW/cm 2 , 7(x) 
is the local evanescent intensity, given by Eq. (2), p(x) 
is the local density, fku is the energy of a photon and 

r = 2tt x 6 MHz is the natural linewidth. The factor X 

15 

appears because linearly polarized light is used. The inte- 
gration is over the entire volume of the evanescent wave. 
Because the absorption is so low, Beer's law remains in its 
linear regime. Obviously saturation effects are taken into 
account. We also account for the Van der Waals potential, 
which leads to a decrease in the atomic density close to the 
surface. Finally also the transverse density distribution of 
the atomic cloud is taken into account. Neglecting these 
last three effects would increase the calculated absorption 
by approximately 20%. The open circles in Fig. 2b are the 
results of this calculation for a density near the prism of 
1.2 x 10 9 cm- 3 . 

For another measurement (not shown) with an angle 
of incidence (9n — c ) = 130 /irad and an evanescent 
probe power of 2.9 /xW, the measured maximum absorp- 
tion of (0.23 ± 0.01)% resulted in a calculated density of 
(1.3±0.4) x 10 9 cm -3 . This value agrees very well with the 
density of (1.3±0.2) x 10 9 cm -3 derived from conventional 
absorption images with a CCD camera. 

Close to a dielectric surface the radiative properties of 
atoms are altered [15]. The natural linewidth of a two- 
level atom can be up to 2.3 times larger than the natural 
linewidth in free space. However, including this effect in 
the calculations only increased the calculated absorption 
by about 2%, which is well within the measurement un- 
certainty. By decreasing both the probe intensity and the 
decay length this effect should start to play a role. 

As a cross check, the second method determines the 
absorption by describing the atomic cloud by a complex 
index of refraction n = 1 + , with an = X4^- the 
resonant cross section for linearly polarized light and p 
the density at the prism surface. Using this index of re- 
fraction to calculate the reflected intensity also yields the 
absorption. The reflected field is determined by evaluating 



E(x,z) 



1 



tt/2 



r p (9, n) exp[ik x (9)x — ik z (9)z 



+inkoUe-0 ) 



2 _ (e-e ) 2 



\d9, (3) 



with k z (9) = nkocos9 the wavevector perpendicular to 
the surface and r p (9) = (cos 9 — n\J 1 — n 2 sin 2 9) / (cos 9 + 
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n\/ 1 — n 2 sin 2 



is the Fresnel coefficient for reflection for 
p polarized light. The same normalization as for Eq. (1) 
was used. The reflected intensity is given by I(x, z) = 
Iq\E{x, z)\ 2 . Saturation effects are not included. Since fi- 
nally only the total absorbed power is important, it is 
not necessary to incorporate the transverse distribution 
in these calculations. The open triangles in Fig. 2b show 
the results of these calculations for various angles of in- 
cidence. The absorption for a maximum density near the 
prism surface of 1.2 x 10 9 cm~ 3 calculated with the com- 
plex index of refraction is slightly higher than the absorp- 
tion calculated from the scattering of evanescent photons, 
mainly because saturation effects were neglected. 



4 Trapping 

In a second experiment we used evanescent waves to de- 
tect atoms trapped close to the surface in a standing light 
field. We load and trap the atoms using the scheme as 
described in previous work [6]. Cold atoms are prepared 
as in the time of flight experiment. During their fall the 
atoms are optically pumped to the F= 1 hyperfine ground 
state. On the vacuum side of the prism surface a repul- 
sive EW potential is created by total internal reflection 
of a 90 mW, TM polarized laser beam with a waist of 
500 fim and blue detuned by 0.2 - 1 GHz with respect 
to the F= 1 <-> F'= 2 transition of the Dl line. This 
potential acts as a mirror for atoms. The decay length of 
the EW can be controlled by changing the angle of inci- 
dence of the laser beam [16]. By scattering EW photons, 
the atoms can make a spontaneous Raman transition to 
the F=2 hyperfine ground state, for which the repulsive 
EW potential is lower. This results in a virtually complete 
loss of their gravitationally acquired kinetic energy [17]. 

The trapping laser is a linearly polarized, 1.3W laser 
beam, red detuned by about 85 GHz with respect to the 
D2 line. It is reflected from the vacuum side of the un- 
coated prism surface, creating a standing wave (SW) with 
a visibility of 0.38. The angle of incidence is nearly nor- 
mal, 25 mrad. The spot size at the prism surface is 380 /im 
x 440 fim (1/e 2 radii). For atoms in the F=l hyperfine 
ground state the EW potential is dominant, whereas for 
the F=2 hyperfine ground state the SW potential is dom- 
inant. Atoms that are Raman transferred to F=2 near 
their motional turning point can be trapped in the SW 
potential. Only a few potential minima will be occupied 
due to the highly localized optical scattering of the EW. 
When the atoms fall towards the prism both the EW and 
the SW are off. Both are turned on 1 ms before the max- 
imum atom density reaches the prism surface. In order to 
decrease the scattering rate, the EW is switched off after 
2 ms, after which the loading of the SW trap stops. 

The EW probe beam is aligned by overlapping it with 
the EW bouncer beam, whose angle can be set with an ac- 
curacy of 25 /xrad. The overlap could be checked over a dis- 
tance of 1.5 m, resulting in an uncertainty in the angle of 
100 /irad. During the trapping experiments the probe an- 
gle of incidence was kept constant at (0q — 6 c ) = 130 /irad 
and the power of the probe beam is Pj n = 2.2 /iW. At 
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Fig. 3. Measurement of atoms trapped in a standing wave, de- 
tected using evanescent probing. The absorption for different 
probe detunings shows that the atoms are distributed over sev- 
eral light shifts. We determine that 1.5 x 10 4 atoms are initially 
trapped. 

the prism surface it had a waist (1/e 2 intensity radius) 
of (770 ± 10) /im. The probe was resonant for atoms in 
the F = 2 hyperfine ground state and was turned on at 
least 14 ms after shutting the EW bouncer beam in order 
to be certain that no falling atoms were present and only 
trapped atoms would be probed. 

Since the trap is red detuned, the atoms will be trapped 
in the intensity maxima. In the center of the trap the scat- 
tering rate in these maxima is calculated to be 7 ms -1 and 
the ground state light shift is <5ls/27t = — 15.4 MHz. 
The trap depth here is only 8.6 MHz since the fringe vis- 
ibility of the standing wave is 0.38. The trap frequency is 
359 kHz, which results in approximately 24 bound levels. 
The resonance frequency of trapped atoms is blue-shifted 
by — 2<5ls, due to the light shift of the excited state. 

Fig. 3 shows the absorption of an evanescent probe 
beam by atoms trapped in the standing wave. The evanes- 
cent wave bouncer is blue detuned by 550 MHz and has 
a decay length £ = 1.15A. The evanescent probing tech- 
nique was used to optimize the trap loading by varying the 
bouncer detuning. The maximum absorption is observed 
for a detuning of the EW probe of 8 MHz. The measured 
linewidth is larger than the 6 MHz atomic linewidth, prob- 
ably due to inhomogeneous broadening. There are two 
contributions to the broadening. Firstly, since the trap 
laser has a Gaussian profile the atoms encounter a spa- 
tially varying light shift. Secondly atoms in higher excited 
vibrational states will encounter a smaller light shift. It 
is not possible to reliably retrieve these detailed distri- 
butions from our measured curve. It is however possible 
to make a good estimate of the total number of trapped 
atoms. The relative absorption curve in Fig. 3 is, for low 
saturation, described by 



-SP 



P 



hwP 
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s(x) p(pc,A) 



1 + 4(4=4) 



d 3 x dA, (4) 



which is similar to the scattering analysis of the falling 
atoms, but inhomogeneously broadened by means of the 
convolution with the Lorentz function. The factor P\„ is 
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the power of the incident probe beam and /5(x, A) is the 
distribution of atoms in the F = 2 hyperfine ground state 
over spatial coordinates and light shifts. By integrating 
Eq. (4) over the detuning S (the area under the curve of 
Fig. 3), the integration of the Lorentzian over 5 yields 
nr/2. The integration over A, for which the integrand is 
now only p(x, A), yields the density of atoms p(x). 

From comparing the kinetic energy of the falling atoms 
in F = 1 to the trap depth when the atoms are pumped 
to F = 2, it follows that mainly the third and fourth 
potential minima will be populated. An underestimate of 
the number of atoms can be obtained by assuming the 
probe intensity constant over the trap region (which is 
valid if the evanescent probe size is much larger than the 
trap size) and all atoms are in the third minimum. Eq. 
(4) then reduces to (7r^r 2 s 3 )/4P in )A^, with N the total 
number of trapped atoms and where S3 denotes the sat- 
uration parameter in the third potential minimum. From 
this a number of trapped atoms in the F = 2 hyperfine 
ground state of 3.0 x 10 3 can be calculated after 14 ms 
of trapping. The total number of trapped atoms will be 
4.5 x 10 3 because the steady state F = 1 groundstate 
population will be 0.5 times the F — 2 population due to 
scattering by the standing wave with the present detuning. 
The populations of the magnetic sublevels of one hyperfine 
groundstate are equal to within 1%. By comparison with 
previous measurements we deduce a lifetime of 14.3 ms. 
An extrapolation results in about 1.2 x 10 4 trapped atoms 
at t = 0. The assumption that the evanescent probe size is 
much larger than the trap size is not completely valid. Tak- 
ing the correct radii into account leads to a 22% increase, 
thus 1.5 x 10 4 trapped atoms. Assuming the transverse 
trap distribution equal to the trapping laser, the vertical 
trap radius to be A/4 and the atoms to be distributed 
mainly over two potential minima, the density becomes 
1.2 x 10 11 cm~ 3 , which is about 100 times higher than the 
density of the atoms falling onto the prism. 

5 Conclusions and outlook 

We have shown that the absorption of a weak, resonant 
evanescent wave can be used to selectively detect cold 
atoms near (~ A) a dielectric surface. A model treating the 
absorption by scattering evanescent photons was suitable 
to describe the absorption. When calculating the evanes- 
cent intensity distribution, the Gaussian character of the 
incident beam had to be taken into account in order to 
quantitatively understand the absorption for angles of in- 
cidence very close to the critical angle. 

By detecting cold atoms falling onto the dielectric sur- 
face for different angles of incidence of the probe beam we 
have verified our quantitative description of the depen- 
dence of the absorption on the angle of incidence of the 
probe beam. By detecting cold atoms trapped in stand- 
ing wave potential minima close to the surface we have 
determined that we have trapped more than 1.5 x 10 4 
atoms initially. This results in an increase of the density 
of two orders of magnitude with respect to the approach- 
ing atoms. 



The technique can be extended to using a CCD cam- 
era so that a transverse distribution of the atoms can be 
measured. By performing measurements for different an- 
gles of incidence of the probe beam, it should be possible 
to obtain information about the distribution of atoms per- 
pendicular to the surface. 

This work is part of the research program of the "Stichting 
voor Fundamenteel Onderzoek van de Materie" (Foundation 
for the Fundamental Research on Matter) and was made pos- 
sible by financial support from the "Nederlandse Organisatie 
voor Wetenschappelijk Onderzoek" (Netherlands Organization 
for the Advancement of Research). R.S. has been financially 
supported by the Royal Netherlands Academy of Arts and Sci- 
ences. 



References 

[I] W. Hansel, J. Reichel, P. Hommelhoff and T.W. 
Hansen, Phys. Rev. Lett. 86, 608 (2001). 

[2] D. Miiller, E.A. Cornell, M. Prevedelli, P. Schwindt, 

A. Zozulya and D. Anderson, Opt. Lett. 25, 1382 

(2000) . 

[3] D. Cassettari, B. Hessmo, R. Folman, T. Maier and 

J. Schmiedmayer, Phys. Rev. Lett. 85, 5483 (2000). 

[4] W. Hansel, P. Hommelhoff, T.W. Hansch and J. Re- 

ichel, Nature 413, 498 (2001). 

[5] H. Ott, J. Fortagh, G. Schlotterbeck, A. Grossmann 

and C. Zimmermann, Phys. Rev. Lett. 87, 230401 

(2001) . 

[6] R. J. C. Spreeuw, D. Voigt, B. T. Wolschrijn and H. 

B. van Linden van den Heuvell, Phys. Rev. A 61, 
053604 (2000). 

[7] Yu. B. Ovchinnikov, I. Manek and R. Grimm, Phys. 

Rev. Lett. 79, 2225 (1997). 
[8] H. Gauck, M. Hartl, D. Schneble, H. Schnitzler, 

T. Pfau and J. Mlynek, Phys. Rev. Lett. 81, 5298 

(1998). 

[9] A. Aspect, R. Kaiser, N. Vansteenkiste, P. Vignolo, 

C. I. Westbrook, Phys. Rev. A 52, 4704 (1995). 

[10] K. Kiersnowski, L. Jozefowski and T. Dohnlik, Phys. 

Rev. A ,57, R4079 (1998). 

[II] K. Zhao, Z. Wu and H.M. Lai, J. Opt. Soc. Am. B 
18, 1904 (2001). 

[12] F. Goos and H. Hanchen, Ann. Phys. (Leipzig) 1, 

333 (1947). 

[13] C. Chiu Chan and T. Tamir, J. Opt. Soc. Am. A 4, 

655 (1987). 

[14] F.I. Baida, D. van Labeke and J.-M. Vigoureux, J. 

Opt. Soc. Am. A 17, 858 (2000). 
[15] J.-Y. Courtois, J.-M. Courty and J.C. Mertz, Phys. 

Rev. A 58, 1862 (1996). 
[16] D. Voigt, B.T. Wolschrijn, N. Bhattacharya, R.J.C. 

Spreeuw and H.B. van Linden van den Heuvell, Phys. 

Rev. A 61, 063412 (2000). 
[17] B.T. Wolschrijn, D. Voigt, R. Janssen, R.A. Cor- 

nelussen, N. Bhattacharya, R.J.C. Spreeuw and H.B. 

van Linden van den Heuvell, Phys. Rev. A 64, 

065403 (2001). 



